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ANALYSIS OF AN INTERFACE CRACK FOR A FUNCTIONALLY GRADED 
STRIP SANDWICHED BETWEEN TWO HOMOGENEOUS LAYERS OF 

FINITE THICKNESS 


N.I. Shbeeb and W.K. Binienda 
University of Akron 
Department of Civil Engineering 
Akron, Ohio 


Abstract - The interface crack problem for a composite layer that consists of a 
homogeneous substrate, coating and a non-homogeneous interface was formulated for 
singular integral equations with Cauchy kernels and integrated using the Lobatto- 
Chebyshev collocation technique. Mixed-mode Stress Intensity Factors and Strain Energy 
Release Rates were calculated. The Stress Intensity Factors were compared for accuracy 
with relevant results previously published. The parametric studies were conducted for the 
various thickness of each layer and for various non-homogeneity ratios. Particular 
application to the Zirconia thermal barrier on steel substrate is demonstrated. 


N AS A/CR— 1999-208874 



1. INTRODUCTION 


Advanced turbine systems and some aerospace applications require use of structural 
ceramics to protect the hot sections. The thermomechanical mismatch between metal and 
ceramics induces high residual stresses responsible for cracking and spallation. One way 
of reduction of the residual stresses is accomplished by processing fully tailored materials 
and interfacial zones with predetermined, continuously varying mechanical properties 
known as Functionally Graded Materials (FGM) [1,2]. FGM could be described as two- 
phase particulate composites where the volume fractions of its constituents differ 
continuously in the thickness direction [3-6], 

Erdogan in his paper [7] discussed the problem of crack growth in FGM due to 
fatigue, creep and stress crack corrosion cracking, and fracture instability. 

He concluded the following: 

1. By eliminating the discontinuity in material property distributions, the 
mathematical anomalies regarding the crack tip stress oscillations for the 
interface cracks are eliminated. Hence, one can now use the crack tip finite 
element modeling developed for the ordinary square-root singularity and apply 
the methods of the energy balance-based theories of the conventional fracture 

mechanics. 

2. Application of FGM as interfacial zones in joining generally incompatible 
materials would greatly improve the bonding strength. 

3. Use of FGM as coatings and interfacial zones would reduce the magnitude of 
the residual and thermal stresses. 
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4 . 


Use of FGM coatings and interfaces would eliminate the stress singularities at 
the points of intersection of interfaces and stress-free ends in bonded materials. 

5. Replacing homogeneous coatings by FGM layers would both enhance the 
bonding strength and reduce the driving forces at the crack tips. 

Delale and Erdogan [8] solved the crack problem for a nonhomogeneous plate. The 
authors considered the plane elasticity problem in which the material is isotropic, has a 
constant Poisson’s ratio (v), and the Young’s modulus (E) is of an exponential form 
varying in the x-direction. They found that the Poisson’s ratio did not have much effect 
on the resulting stress intensity factors. They also found that the strain-energy release rate 
of the crack embedded in the portion of the medium with higher stiffness is lower than 
that corresponding to the crack tip in the less stiff side of the material. Hence, the crack 
will grow in the direction of the less stiff material. 

Delale and Erdogan [9] considered the interface crack in a nonhomogeneous elastic 
medium. In this paper the interface crack between two bonded half planes was addressed. 
One of the half planes was homogeneous while the other was nonhomogeneous in a 
manner that the elastic properties are continuous throughout the plane and have 
discontinuous derivatives along the interface. They assumed that the Young’s modulus 
and the Poisson’s ratio are of the exponential form. They found that the singular behavior 
of the stress state near the crack tip in the nonhomogeneous medium is identical to that in 
a homogeneous material given that the spatial distribution of the material properties are 
continuous near and at the crack tip. 

Also, Delale and Erdogan [10] solved the crack problem for two bonded 
dissimilar homogeneous elastic half-planes and assumed that the interfacial region, can 
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be modeled by a very thin non-homogeneous layer. The elastic properties of the 
interfacial material varied continuously between that of the two semi-infinite planes. It 
was assumed again that E and v are exponentially dependent on y. Varying the Poisson’s 
ratio did not have much effect on the stress intensity factors. The Airy stress function was 
used in their formulation of the solution in which it was assumed that it is composed of 
two functions, one is associated with an infinite plane containing the crack on the x-axis, 
while the second is an uncracked strip. Their results showed that if the crack location 
approaches the less stiff material the strain energy release rate increases. 

Erdogan et al. [11] studied the perpendicular crack to the interface in a bonded non- 
homogeneous material. The main goal was to study how the singular behavior of the 
stresses and stress intensity factors is affected by the very steep variations in the material 
properties near the diffusion plane. Here they assumed that the shear modulus (p) varies 
in an exponential form. They suggested that: 

1 . Regardless of the mechanism of binding at the atomic level, in many cases there 
is always a thermodynamically stable and readily distinguishable region 
between the two homogeneous materials. 

2. The interfacial regions are generally locations of higher concentration of stress 
and micro-flaws. 

3. In most material pairs the fracture toughness and the sub-critical crack growth 
resistance of the interfacial zone tend to be lower than that of adjacent 
homogenous materials. 

They found that the nonhomogeneity constant P has a great effect on the stress 
intensity factors. They showed that as P increases so do the stress intensity factors. 
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Erdogan and Ozturk [12] solved the mixed boundary value problem for a non- 
homogeneous medium bonded to a rigid subspace. They investigated a two-dimensional 
diffusion problem in which the interface contains a plane crack. Rather than solving the 
problem for a given material, they used an inverse method. This inverse method provided 
the material constitutive behavior for which the mixed boundary value problem could be 
solved. 

One year later, Erdogan and Ozturk [13] studied the axisymmetric crack problem in 
a non-homogeneous medium. They noticed that the crack opening displacement in non- 
homogeneous materials was significantly greater than the corresponding homogeneous 
values. 

Later Konda and Erdogan [14] considered the mixed mode crack problem in a non- 
homogeneous elastic medium. The crack was arbitrarily oriented with respect to the 
direction of the property gradient. The effect of the variation of the Poisson’s ratio was 
neglected in the solution because of previous studies. 

Erdogan and Wu [15] studied the crack problem in FGM layers under thermal 
stresses. They considered an unconstrained elastic layer under statically self-equilibrating 
thermal or residual stresses. The layer contained an embedded or surface crack 
perpendicular to its boundaries. After giving the distribution of thermal stresses, the stress 
intensity factors for the embedded and surface crack were presented along with the 
results of the crack/contact problem in a FGM layer that was under compression near and 
at the surface and tension in the interior region. 

Chen and Erdogan [16] solved the interface crack problem for a nonhomogeneous 
ceramics coating bonded to a homogenous metallic substrate using displacement 
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formulation. It was concluded in this work that the dominant mode of the stress intensity 
factor is Mode I, and it decreases as the nonhomogeneity constant changes from soft to 
stiff under uniform normal stress. A similar trend was noticed for Mode II, under uniform 
shear stress. Decreasing the thickness of either material increased the stress intensity 
factors. 

In this work the problem described in Chen and Erdogan [16] will be extended to 
include the third thin layer of homogeneous ceramics material to increase thermal 
protection of the metallic substrate. As shown by Kokini and Choule [17], thermal 
barriers always include some thickness of pure ceramics material. The scope of this work 
includes an examination of debonding of ceramics layers from the substrate. In particular, 
the stress intensity factors (SIF) and strain energy release rates (SERR) are obtained 
using the Airy stress function formulation for the interface crack embedded between the 
finite thickness substrate and the non-homogeneous strip, which is sandwiched between 
the substrate layer and thin homogeneous layer. It is assumed that the FGM has a 
constant Poisson’s ratio and the shear modulus is of an exponential form. Plane elasticity 
is assumed and the solution is valid for both plane stress and plane strain. 
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2. FORMULATION 


The geometry of the problem is shown in Figure 1. Both the substrate and the 
coating, which are perfectly bonded to the FGM, are isotropic and homogeneous, and 
have h| and h} as their respective thickness. The FGM thickness is hi, and is denoted as 
material 2. 

In the global x-y coordinates the shear modulus of the FGM is assumed to be as 
follows: 

H 2 (y) = H ] e r> (1) 


where 


7 = 

h, p. 


( 2 ) 


Hooke’s law relates strain and stress using two independent material constants: 


ejx, y) = = -M(k + iKv +(^-3Kvl 

dx 8 // 

£ n (x, y ) = |^- = -L [(fc - 3>y„ +{k + 1>t vv ] 
ay 8 p 

1 du dv 1 


(3) 


where u and v are horizontal and vertical displacements, and k is defined as: 

k = 3 - 4v for plane strain 

K = - — — for plane stress 

1+v 
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The Airy stress function F(x,y) is defined by 


cx„(.*,y) = 
o - vv U,y) = 
T (V U,y) = - 


d 2 F(x,y ) 

d^Fj^y) 

dx 2 

d^F\x, y) 
dxd y 


(4) 


By incorporating (4) and the strain-stress relations in the following compatibility 
equation: 


d 2 e d 2 £ d 2 e 

— ~^ + f-2 — = 0 

d .V d y~ dxdy 

we obtain the fourth order governing equation for F(x,y): 

d 4 F ] {x,y) | ^ <9 4 F,(x, y) | d^j^y) Q 


dx 4 dx 2 dy 2 


d v 4 


The solution of (6) is found by applying the Fourier Transform: 


V(a,y)=]F i (x,y)e ,a 'dx 


(5) 


( 6 ) 


(7) 


that transforms (6) into 4-th order Ordinary Differential Equation (ODE) with 
constant coefficients. The ODE is solved by forming the characteristic equation to 
determine its roots: 


4l/ . ,d 2 V d 4 V n 
a V -2a- — r + — - = 0 
dy 2 dy 4 

=> w 4 -la, 2 w 2 + a 4 =0 

=> w, = h', = |a| & w 2 = w 4 = -|a| 


( 8 ) 
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Thus, the stress function for the substrate has a form with double repeated roots as follows: 


t 00 

F x ( x , y) = — f [(£>, (a) + yD , (a))e |flr|v + (D,(a) + yD 4 (a) )e >|v ] e ,xtt da (9) 
2k j 


Similarly, for the coating. 


F i (x,y) = —][(C l (a)+yCAa))e laly +(C,(a)+yC 4 (a))e' la]y }e ixa da (10) 
2k j 


The governing equation for the FGM is: 


V 2 ((7 „+cr vv ) + r 

V 4 F 2 (x,y) + y 2 


K-, -3 

CT » + ^TT a ' 

JC-, + 1 


- 2 Y-K-«r» + CT vv) = ° 

^ v 


d 2 F 2 (x,y) + k 2 - 3 d 2 F 2 (x, y) 


d y 2 fc 2 +l d x 2 


(ID 


-2y— V 2 F,(x,>9 = 0 

ay 


The characteristic equation of ( 1 1 ) becomes, 


m 4 -2ynr +(/ 2 -2 a 2 )m 2 + 2 ya 2 m + (a 4 -a 2 y 2 — — -) = 0 

K" 1 I 1 


( 12 ) 


There are four independent roots of (12): 


3- k 2 
k 2 +1 

3-k 2 

K i + 1 

3-k 2 

K, + 1 

3-k 2 

2 ' V** ' 4 ‘^V^+l 


w, = T-,l«‘ +i r+*a)\ 




a" +— +/ay, 

4 M 


w, =t -Ja‘+- L r- , «y, 


7 


rw — — — -1- TV " -I- / 7YV 


Consequently, the stress function is obtained in the following form: 


FAx, v) = — [4, (a)e'" 1 ' + A,(a)e" ! - V + AAoc)e'"" + A 4 (a)e" ,jV ] e' to Ja 
- - 2k j 


(13) 


(14) 
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( 15 ) 


Using (4) the stresses for the FGM are: 

cr^ix, v) = — f[m, 2 V m,v +m;A,e m >' +m;A,e"" x + m 2 4 A 4 e n '* y ]e ixa da 
2 n 

a [; ) (x,y) = -— \a 2 [A l e’"' y +A 1 e m '- y +A,e"‘ x ' + A 4 e m " ]e Ua da 
2k 

r' 2) (x, y) = — — 1 ia[m ] A l e m '' +m- l A 1 e"’ 2 ' +m 4 A 4 e nh ' + m 4 A 4 e" u ' \ e‘ m da 
2k j " 


the stresses for the substrate are: 

I 00 

tOx, y) = — J[(or(£>, + yD 2 ) + 2\a\D 2 )e H ' +(cc 2 (D , + yD 4 )-2\a\D 4 )e My ]e im da 


<‘U.V) = ~ \a 2 [(D t + yD ,)] e |a|v +(D, + yD 4 )] 
2/r J 


(16) 


C (.v, y) = - 2 - | /a[(|a|(D, + yD 2 ) + D 2 )e My + (-|a|(D, + yD 4 ) + D 4 )e M ' ]e ,xa da 


and for the coating they are as follows: 

< ( x, y ) = 2- j[(or (C, + yC 2 ) + 2\a\C,_ )e W ' + (a 2 (C, + yC 4 ) - 2\a\C 4 )e M ' }e ixa da 


2k 

\ 

2k 


<CU >’) = -2- (a 2 [(C, + yC 2 )k |a|v + (C, + yC 4 )e~ Hx }e lxa da 

/ TT J 


(17) 


T^{x,y) = 


_1_ 

2;r 


J /a[(|a|(C, + yC 2 ) + C 2 )e ]al ' + (-|a|(C 3 + yC 4 ) + C 4 )e-‘ |a| v ]<?“"</« 


From (15), (16) and (17), it can be seen that there are 12 unknown constants (in the 
Fourier space they are functions of a), i.e., D s , Q and Aj (i=1...4), which can be obtained 
from the following conditions: 

O^(x-h l ) = 0,t' x [ ) (x-f h ) = 0 
o i2) {x,0) = a;' v ’(.v,0), t ( , 2 ’(x, 0) = C(x,0) 
a £’(*, h 2 + h 2 ) = 0, (x, h 2 + h 2 ) = 0 

o {2) (x, h 2 ) = <r'f(x, h 2 ), t‘;’(x, /i 2 ) = (x, h 2 ) 

u ( 2) (x, h 2 ) = u°'(x, h 2 ), v t2) (x, h 2 ) = v ,?) (x, h 2 ) 

- °0 < * < OO 


( 18 ) 
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and the following mixed boundary conditions: 


a <2) (x,0 + ) = p, (x), t ( 2) (x,0 + ) = p 2 (x), -a < x < a (19a,b) 


u' 2) (x,0 + )-u <l, (x,0 ) = 0, v <2, (x,0 + )- v"'(x,0 ) = 0 |xj > a (20a,b) 

where pi and P 2 are known surface tractions. In order to convert the solution from a dual 
integral path into a singular integral equations, (20a,b) has to be replaced by the 
following auxiliary functions: 


— [u <2, (x,0 + ) — u (1) (x,0 )] = f,(x) -oo< x <o° 

dx 

— [v l2 ’(x,0 + ) - v (l, (x,0“)] = f,(x) -oo < x < °° 

3x 


(21a,b) 


Consequently, all twelve unknowns are determined in terms of the auxiliary 
functions, f|(x) and f 2 (x), by using (18) and (21). The unknown auxiliary functions are 
solved by using (19). After some lengthy manipulations the following singular integral 
equation (SIE) was reached. For details see [18], 


7r(l + x,)(l + K: 2 ) 


^(-A(Jt))= ]^dt + \f,{t)K u {xj)+\f z (t)K v \x,t) 

i “f - /C* t ) t X 


(2 + K*j + K* 2 ) 


[ 7ry(5 4- 5fC 2 -f ^ + *7 (2 + ?c 2 ) + k x (6 + 4k 2 )) Ut-x) 

4(1 + tCj )(1 + jc 2 )(2 + x - ! + x 2 ) i \t-x\ 2 

- y (1 ZJSlZJSl 1 2k > ( 2 + k 2 ) + K{(2 + kJ) f Ci(U (t - x))j\{t)dt 
2 ( 1 + x , )(1 + k - 2 )(2 + x , + k 2 ) 4 


7T(l + K-,)(l + y 2 ) ( _ ft (JC)) = j A +jf } {t) K 2x (X, /) + } f 2 (i t)K 22 ( A- , t) 
/I x (2 + K l + K 2 ) -a t ~ X -a -u 

| ^■y(l + 2y 2 + K: 2 +x- 2 2 (2 + y 2 ) + x- | :! (2 + x- 2 )) jr (/-j:) 

2( 1 + x, )(1 + x 2 )(2 + x, + x 2 ) 

y( 1 - x, - k ; + 2X|(2 + x 2 ) + x 2 (2 + x 2 )) 


( 22 ) 


+ 


2(l + x,)(l + x 1 )(2+x l + x 2 ) 


jCi(U{t-x))f 2 (t)dt 
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where. 


AT n (x,/) = 


(1 + K - , )(1 + K" 2 ) 


2(2 + K - , + k 2 ) o 


{](•*,,« - x u )ism(a(t-x))da + 


t' 

J(x n + .r ll( -2( 


(1 + JC, )(1 H- #c* 2 ) 


)) cos(a(/ - x))da + J (jc, , 


+ x 


lie 


(23) 


- 2 ( 


V, — K ^ 


(1 + K) )(1 + K\ ) 


, y(l + 4jf, +2k s -k\ +2v,v, + , 

+ r— 777 : ="7 s -))cos(a(r - x))da } 

2g'(1 + k’ 1 )~(1 + k-,)~ 


K,,(x,t) = 


(1 + k,)(1 + k 2 ) 
2(2 + K, + K, ) 


{J(X|2c -X 12 +2i( 


2+ K, +k 2 
(1 + k,)(1 + k 2 ) 


y(5 + 6k, + 2k, 2 + 5k, +4k,k, + k, 2 k, + k; . . 

1 1 — — — -) i sin(a(t - x))da 

2oc(1 + k,)' 1 (1 + k 2 ) 2 

+ J(x l2 + x 12c )cos(a(t-x))da} 


(24) 


K„(x,t) = 


(1 + K, )(1 + K 2 ) 

2(2 + k, +k 2 ) 'l 


( J( x 21c — X 21 +2i ( 


2+ K, + K, 

(1 + K, )(1 + K, ) 


Y(1 + 2k,+2k, + k, + k,k, + k: . . 

- 1 - — L — t— 1 77" ■-) 1 sin(a(t - x))da 

2a(l + K,)-(l + K,)- 

+ J(x 2l + x 2k .)cos(a(t-x))da| 


(25) 


K^(x,t) = 


(1 + K, )(1 + K, ) 


2(2+ k, + k 2 ) J 


{J(x 22c -x 22 )isin(a(t-x))da + 


J(x 22 + x 22c -2( 


k, — K\ 


(1 + K, )(1 + K,) 


))cos(a(t - x))da + J (x, , 


+ X 


lie 


— 2 (- 


K, — k 2 Y(1 + 4k, + 2k 2 -k 2 +2k,k 2 + k, 2 k 2 - k 2 


(26) 


(1 + K, )(1 + K, ) 


w, I — IX] ivt I I IX, IV v , / /* UJ 1 

. V, -))cos(a(t -x))da} 

2a(l + K,)'(l + K,)' 
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C 0 is the Euler constant, whereas U is the upper limit beyond which the cosine 
integral, Ci, is negligible. Q is the determinant of the 4 by 4 coefficient matrix and Qy are 
the corresponding 3 by 3 cofactors. Both can be found in the Appendix. 
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3. SINGULAR INTEGRAL EQUATIONS SOLUTION 

The singular integral equations (22) contain Cauchy kernels. In order to obtain 
unique results the following conditions need to be incorporated into the solution. 


a 

jf,U)dt = 0 i = l,2 


(29) 


The system of equations (22) and (29) can be solved together using the Gauss 
quadrature method. For example, using, the Lobatto-Chebyshev collocation method was 
shown in [19] to produce accurate results for the above type of equations. In practice the 
auxiliary functions are discretized at particular points t k , integration is replaced by a 
summation, and the system of linear algebraic equations are obtained for collocation 
points x p (x p is never the same as t k ) in the following form: 


h g . ( tj. ) vi^ 


>=i it *=i t k x p j= 1 1= , 


(30) 


where p=l,...,n, w k is the weighting coefficient, and R„ is a remainder that becomes small 
for sufficiently large number of points t k . According to the aforementioned technique, the 
abscissas are calculated according to: 


,(*- 1)jt x 

t k = COS(- ; — ) • 


n — 1 


•k = 1 


(31) 


The corresponding weights are: 


w i = w „ = 


k n 

;h' = 

2(n — 1) n — 1 


• r = 2, . . . , n — 1 . 


(32) 


The collocation points can be found from: 


X p = cos( p- — — ) p = \,...,n-l. 

2 / 1-2 


(33) 
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The remaining two equations are generated using (29) in the following form: 


(s k )w t = 0 

Jt=l 

n 

^8i( s k) w t = ° 


( 34 ) 


k=\ 


By incorporating (30), (34) and (33) together, the system of equations can be represented 
as follows: 

[AL, ! „feL={f’L (35) 


thus, the unknowns are obtained by: 




(36) 


The mode-I and mode-II SIF are defined as follows: 


k,(a) = lim > /2(x 1 -a)o yiyi (x,,0) 


(37) 


k 2 (a)= lim^2(x, -a)T X|) . (x,,0) 

From the principal part of the expressions for gi(t) and g 2 (t) as shown in [20], the 
following are obtained for ki(a) and k 2 (a): 

1 -i_ V ^ 

(38) 


|X| (2 + K, + K-, ) 

k,(a) = -^ — — 1 ^7=g 2 ( a ) 


(1 + k, )(l + K,)Va 


, , , JI,(2+ K, + K,) 

k,(a) = — — — — 1 ^=g,(a) 


(39) 


(l + K,)(l + K 2 )Va 

Where g,(a) and g 2 (a) correspond to gi(l) and g 2 ( 1) respectively when solving (36). 
The strain energy release rates (SERR) can be calculated from [16], and they are 
listed as follows: 


G,(«) = — ^ A'i(a) 


G-,(a) = 


8/t, 

jz{k 2 + 1 ) 

8jU, 


(40) 
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where G ( and G 2 are the opening mode and sliding mode SERR, respectively. The total 
SERR is expressed as: 

G T (a) = — (a) 2 + k-, (a) 2 ) (41) 

8//, 

Before performing the parametric study the above solution should be verified with an 
existing solution close to the proposed problem. The verification is accomplished by 
comparing the results of our model with that of [16] by letting h 2 approach zero and ya 
(the normalized non-homogeneity constant with respect to the crack length) taking the 
values as shown in the Table I. 

The values were obtained for hi=2a, h 2 =a, V|= v 2 = v 2 = 0.3 under loading of uniform 
normal stress. The results in Table 1 demonstrate the accuracy of the solution. 

4. RESULTS 

The geometry of the problem being examined is shown in Figure 1. The thickness 
of each layer is normalized with respect to the half-length of the crack “a” located at the 
interface between the substrate layer and FGM. The homogeneous substrate, material 
“1”, may be stiffer or softer with respect to the homogeneous layer of ceramics, material 
“3”. The normalized nonhomogeneity constant ya is varied between -3 and 3, which 
covers most of the practical cases. The results are calculated for normalized mode-I and 
mode-II SIFs, i.e., k,/ko and k 2 /k<), and normalized SERR, i.e., G|/G 0 and G 2 /G () , where 

k 0 =<7 Jc and G 0 = — . The results are calculated for two loading conditions, 

7T(K\ + 1) 

namely unit normal stress in y direction and unit shear stress in x-y plane. The Poisson's 
ratio is assumed to be v = 0.3 for each layer. 
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First, the influence of the additional ceramics layer will be examined in order to 
determine if there is any advantage to fabricate a thermal barrier with pure ceramics on 
top of the FGM layer. Figures 2 and 3 show normalized total SERR versus the thickness 
of the homogeneous ceramic layer h ? produced by normal stress and shear stress, 
respectively. The thickness of the substrate is assumed 4 times higher than the thickness 
of the FGM. The crack length is assumed to be the same as the thickness of FGM layer. 

It can be noticed that SERR is significantly reduced even by a small additional 
thickness of the ceramics material. The rate of reduction is the highest for hj below 0.5a 
for h} larger than 3a, the influence is negligible. The largest reduction is produced for the 
negative nonhomogeneity constant when the stiffness of the ceramics material is smaller 
than the substrate. The smallest reduction is in the case of the stiffer ceramic material. 
The homogeneous case is obtained for ya = 0 and its SERR curve is located in between 
the other two cases shown. One can conclude that the small change of the thickness of the 
layer above the crack can significantly reduce the SERR. It can be recommended that the 
optimum thickness of the ceramics layer should be about half of the thickness of the 
FGM layer. 

Using the above recommendation we can compare the SERR with the results 
obtained by Chen and Erdogan (1996) indicated in the following figures by h} = 0 (this 
case is also the limiting case of our model). Figures 4, 5 and 6 are generated for the case 
of the normal applied stress for three different thickness of the FGM layer. The thickness 
of the substrate is assumed to be very large. Comparison of the curves with the additional 
layer of ceramics and without that layer (h 3 = 0) shows that the ceramics layer reduces the 
SERR for each case, but the most significant reduction is for the thinnest FGM layer. In 
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all cases, the magnitude of G, k| and k 2 decreases as the stiffness of the ceramics layer 
increases. 

Figures 7-9 show SERR and SIF for the case of the applied shear stress for the 
same geometrical and material conditions as in the Figures 4-6. Similarly as for normal 
applied stress, G, ki and k 2 decrease with the increasing nonhomogeneity constant. For 
the case of the softest ceramic layer, additional thickness h 2 = 0.125a of the homogeneous 
layer added to 0.25a thick FGM reduces G, k| and k 2 as much as increasing the thickness 
from 0.25a to 0.5a of the FGM. Hence, either increasing the FGM thickness or increasing 
the ceramic layer can reduce the SIF and SERR. 

Figures 10 and 11 show the normalized total SERR generated by normal and 
shear stresses, respectively, versus the nonhomogeneity constant for the various thickness 
of the substrate. The reduction of the SERR by increasing thickness of the ceramic layer 
is equally significant to each thickness of the substrate. The cases for hi = 10a and h] = 4a 
overlap for the applied normal stress and they are almost identical for the applied shear 
stress. As in the previous cases, the stiffness ratio of the ceramics coating to substrate 
significantly changes the SERR, especially for the thin layer of the substrate. 

Mode-I is dominant under normal stress tractions. A question arises what is 
better under such loading conditions: thicker FGM layer without any homogeneous layer 
of ceramics or thinner FGM with h 2 making the difference. In the last parametric study 
we will assume that the total thickness of the FGM and homogeneous ceramic coating is 
constant. The results of SERR are plotted with respect to ln(pV|j.|) in Figure 12. 

It can be noticed that for increasing In(pVpi) the SERR decreases. The highest G 
is for the most negative ln(p 2 /pi). By replacing part of FGM by pure ceramic material we 
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increase G for ln(p.3/|i|)=-2 (the region above the crack is softened by the softer ceramics 
layer) and decrease G slightly for ln(|i3/|J,i)=2 (the region above the crack is stiffened by 
the stiffer ceramics layer). 

Hence, it can be concluded that SERR can be reduced by several methods: stiffer 
coating application, thicker FGM layer, and additional layer of homogeneous ceramics. 
The most optimum combination depends on the stiffness ratio of the ceramics with 
respect to the substrate. 

Let’s apply the above knowledge to the specific cases shown in Figure 13 (also 
discussed in [17]). The substrate material is steel. Zirconia is used for the ceramic 
coating. The FGM is made by gradual change from 100 % of Zirconia to 
Zirconia/CoCrAlY to the bond coat attached to the substrate. In case (a), there are 50/50 
of Zirconia/CoCrAlY layer and bond layer sandwiched between pure ceramics and 
substrate material. In case (b), there are four layers of Zirconia/CoCrAlY, i.e., 75/25, 
50/50, 25/75, and bond layer sandwiched between pure ceramics and substrate material. 

The properties of materials are taken from [21]. For Zirconia they are: E 3 = 36 
GPa and v 3 = 0.2, while for steel they are: E, = 207 GPa and v, = 0.33. The thickness of 
each layer can be normalized with respect to crack length, which is assumed to be the 
same as half of the total thermal barrier thickness, i.e., h, = h 2 +h 3 = a. The thickness of 
the substrate is hi = 6a. 

Table II shows results for three different cases under normal stress conditions. In 
the first case the ceramic layer is part of the FGM. In the second and third cases the pure 
ceramic layer belongs to material “3” and FGM thickness is measured from the interfaces 
between steel and bond coat to the interface between Zirconia and FGM layer. It can be 
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noticed that the total SERR is higher in the lower two cases than it is in the first case by 
16% and 38%, respectively. Hence, it is critical how FGM is defined in the practical case 
studies. It can be also noticed that the thinner Zirconia layer produced smaller SERR 
because of the negative nonhomogeneity ratio as described in Figure 12. Finally, the 
results show some sensitivity to the FGM Poisson’s ratio, but for all practical purposes it 
can be assumed to be the same as for the pure ceramics. 
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Table I 

Verification of the model 


ya 

ki(a)/(G 0 Va) 
Chen and 
Erdogan 
(1996) 

k|(a)/(o<Wa) 

Present 

study 

k 2 (a)/(GoVa) 
Chen and 
Erdogan 
(1996) 

k 2 (a)/(a 0 Va) 

Present 

study 

-3.0 

2.430 

2.428 

-0.681 

-0.624 

-2.5 

2.252 

2.251 

-0.571 

-0.533 

-2.0 

2.087 

2.087 

-0.471 

-0.445 

-1.5 

1.936 

1.936 

-0.379 

-0.364 

-1.25 

1.866 

1.866 

-0.336 

-0.325 

-1.0 

1.799 

1.799 

-0.296 

-0.288 

-0.75 

1.735 

1.735 

-0.258 

-0.252 

-0.5 

1.675 

1.675 

-0.221 

-0.218 

-0.25 

1.618 

1.618 

-0.187 

-0.186 

-0.01 

1.566 

1.566 

-0.156 

-0.155 

0.25 

1.514 

1.514 

-0.125 

-0.126 

0.5 

1.466 

1.466 

-0.096 

-0.099 

0.75 

1.422 

1.422 

-0.069 


1.0 

1.380 

1.379 



1.5 

1 .304 

1 .303 

0.002 

0.003 

2.0 

1.237 

1.237 

0.042 

0.036 

2.5 

1.179 

1.178 

0.077 

0.070 

3.0 

1.128 

1.127 

0.108 

0.100 
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Table II 

SIF and SERR for Zirconia thermal barrier and steel substrate. 


Geometry and 

Material 

Property 

k,(a) 

k 2 (a) 

Gj(a) 

V 2 =V3 

V 2 =V| 

v 2 =v 3 

V 2 =Vi 

v 2 =v 3 

V 2 =V| 

h,=a, h 2 =a, h 3 =0 
ln(|i 2 /|i|)=- 1.646 
ya=- 1 .646 

1.9471 

1.9178 

0.4102 

0.2888 

3.9598 

3.7613 

h,=a, h 2 =0.771a, 
h 3 =0.229a 
ln(ji3/pi)= -1.646 
Ta=-2.135 

2.0905 

2.0336 

.4827 

.3190 

4.6033 

4.2371 

h,=a, h 2 =0.545a, 
h 3 =0.455a 
ln(p 3 /|ii)= -1.646 
ya=-3.021 

2.243 

2.1690 

.5565 

.3448 

5.3389 

4.8236 



Figure 1. Geometry of the interface crack for a functionally graded layer sandwiched 
between the homogeneous substrate and coating materials. 
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Figure 2. Influence of hVa on the total SERR for hi/a=4.0 and h2/a=1.0 
under loading of uniform normal stress. 

1.8 -r-r- n T T — 1 ^ T * T - 1 T 


1.6 



*• — ya=-3.0 
-■ — ya=-0.01 


ya=3.0 



0.6 — 1 1 “ — 1 ‘ *■ 

0 0.5 1 1.5 2 2.5 3 3.5 4 

h/a 

Figure 3. Influence of h-^/a on the total SERR for h|/a=4.0 and h2/a=l .0 
under loading of uniform shear stress. 
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Figure 4. Influence of IWa and hVa on the total normalized SERR for h|/a= 100.0 
under loading of uniform normal stress. 



ya 

Figure 5. Influence of lWa and hVa on the normalized mode I SIF for h|/a= 100.0 
under loading of uniform normal stress. 


N AS A/CR— 1999-208874 


26 




under loading of uniform normal stress. 



- 3 - 2-10123 

ya 

Figure 7. Influence of IWa and hj/a on the total normalized SERR for h|/a= 100.0 
under loading of uniform shear stress. 
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Figure 8. Influence of h 2 /a and h 3 /a on the normalized mode I SIF for hj/a= 100.0 
under loading of uniform shear stress. 
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Figure 9. Influence of h 2 /a and h 3 /a on the normalized mode II SIF for hi/a=100.0 


under loading of uniform shear stress. 
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- 3 - 2-10123 

■ya 

Figure 10. Influence of h]/a and h 2 /a on the total normalized SERR for h 2 /a=l .0 
under loading of uniform normal stress. 



ya 

Figure 1 1. Influence of hi/a and hVa on the total normalized SERR for h 2 /a=1.0 
under loading of uniform shear stress. 
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Figure 12. SERR versus ln(|a. 3 /jJ.i ) for l^+h-* = 0.25a 
under loading of uniform normal stress. 
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h.i 


h, 


h 2 


h. 


Zirconia 


50/50 


Bond Coat CoCrAlY 


Steel Substrate 


0.1018 


0.1018 


0.02 cm 


1 .27 cm 


(a) 


h, 


Zirconia 



75/25 



50/25 

h 2 


25/75 



Bond Coat CoCrAlY 

h, 


Steel Substrate 


0.0509 cm 

0.0509 cm 
0.0509 cm 
0.0509 cm 

0.0191 cm 


1 .27 cm 


(b) 

Figure 13. Geometry of the Zirconia thermal barrier on steel substrate (from Kokini and 
Choules (1995)), (a) two layer FGM, (b) four layer FGM. 
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APPENDIX 


CU< (a ((l-E 4hl<, *4E J1<l ‘ , 'iaa) ♦ 

a (l*i thlfl -2E 2hl0 (l + 4hl 2 a 2 ) - (-1* E 2hla ) 2 K 1 ) } ) / 
(-4 (-1+ E 21 * 3 ) 2 * 16E? W,a hi 2 a 2 ) ♦ 

(-a 2 (-3 +x 2 ) + n£ <l+x 2 )); 

o 

02 » {a ((l-E 4hl0 + 4E? hla hla) nb(l ♦ *l) + 

a (l*g* hl<, -2S? lao (l+4hl 2 a 2 ) - (-1 + E 2hla ) 2 x 1 ) ) ) / 
(-4 (-1+ E Jhla ) 2 + 16 E 21110 hi 2 a 2 } + 

7 (-« 2 (- 3 + * 2 ) ♦ ( 3 + K 2 > ) 1 

o 

03- (at ({1- E 11 * 0 * 4E? hla hla) i % <l + xi) -* 

a (1+ E 41 * 1 ” - 21? 1110 (1 + 4 hi 2 a 2 ) - {-l*E IBlo ) 2 x 1 ) ) ) / 
(-4 (-1+ E? ba ° l ) 2 ♦ 16 E 21110 hi 2 a 2 ) * 

-7 {-“ 2 (-3-I- x 2 ) +xn| (l + x 3 )) ; 

0.4= (a ((l-E IKla + 4E ,hlo hla) m* (1+xi) + 

a (l + fi llac, -2E 2hla (1 + 4hl 2 a 2 ] - (-1+ E 2Hla ) 2 x 1 ) ) ) / 
(-4(-l + E m °) 2 4l6l? hlo hl 2 a 2 ) + 

— (-a 2 (-3 + xj) ♦Utf (l+x 2 )) ; 

8 

C21= -(la (-a (-1+ E lhl “ + 4E ihl “hla) (l + xj.) + 

ll4(l+S lhl “-2E 2hla (1+ 4 hi 2 a 2 ) - (-1+E ? Kla ) 2 x 1 )))/ 

(-4 (-1+ E 2 * 10 ) 2 * 16 8?“" hi 2 a 2 ) + 
la (mi? (-3»x 2 ) - a 2 <l»x 2 )) 

8 (-r -uni) 

C22 - - (la {-a (- 1* E 41 * 0 + 4 E? hla hla) (1+xi) ♦ 

HU (l* E* hlo -2E mo (1+ 4hl 2 a 2 ) - (-l + E Jhl “) 2 x 1 )))/ 
(-4 (-1* E 21 * 0 ) 2 * leE^h^a 2 } + 

Xa (m? (-3 ♦ x^) - a 2 <l+x 3 )) __ 

8(-r+«b) 

C23 • -(la (-a (-1+ E 41 * 10 ♦4B? hla hla) (l + x!> ♦ 

life (!♦ & hLa - 2B? hl0 (1* 4hl 2 a 2 ) - (-1* E 2M “) 2 x 1 )))/ 
(-4 (-I* E 2 * 10 ) 2 * 16 E 21 * 0 hi 2 a 2 ) + 

Xa (m| (-3 ♦ x 2 ) - a 2 (l+x 2 )) 

8 <-r *nfe) 

C24- -(la (-a (-1+ E 41110 ♦ 4 E 21 * 0 hi a) (l*x r ) ♦ 

m t (l*E 11,la -2E 2hla (1+4 hi 2 a 2 ) - (-1* E lhla ) 2 xi) ) ) / 
(-4 (-1+ E 21 *") 2 * 16 E 22 * 01 hi 2 a 2 ) ♦ 

Ia (gj (-3* x 2 ) - a 2 (1» x 2 ) ) _ 

8(-r+m») 
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C31- (IE? a “ 1 ' { (1- E lh3 “ + 4E?“ a h3o) rai (l+*3) ♦ 

a (l*E* Ma -2E?“- (l.ihfa*) - (-1 + E? M “) 2 k 3 )))/ 
(4((-l*E 2 “ a ) 2 -4E?“ a h3 2 a 2 )) ♦ 

I E? 12 **- (a 3 (-3 + X 2 ) - »? (1+ x 2 ) ) . 

8a 

032 = (IE 13 " 2 ( (1- EJ 410 ® + 4E? K3a h3a) nt <1+ x 3 ) + 

a (1+4 h3 2 a 2 ) - (-1+ E?““) 2 k 3 ) ) ) / 

(4{{-l+E ,h3a ) 3 -4E 2 “ a h3 2 a 3 )) ♦ 

( a 2 (_3 +x 2) - ( l+*3>) . 

, . # 

8a 

C33 r (iE“"3((l-^“°+4E 2h3 “h3a) Ufe(l+X3) ♦ 

a (l+fi 4h3a -2E? h3a (l+4h3 2 a 2 ) - (- 1+ E 2 *"}^) ) ) / 
(4((-l+E 2 “ a ) 2 -4J*“ a h3V)) + 

IE h2 “3 (a 2 (-3 + K2) - I"? (!+«)) . 

8a 

C34 - (IE? 12 ** {{1- E?“ a + 4E 2 “°h3a) n>i<l + x 3 ) + 

a (l+E* Ma -2B?“ a (1 + 4h3 2 a 2 ) - (-1+ E?“*) 2 x 3 ) ) ) / 

(4 ((-1+ E? h3o ) 2 -4E 2Mo h3 2 a 2 }) + 

(a 2 (-3 + * 2 ) (l+*3)) 

8a 

041= (E* 121 ^ (-a (-1+ S 4h3a + 4E? h3a K3a) (l+x 3 ) ♦ 

n* (1+ 21?““ (l+4h3 2 a 2 ) - (-1+ I?* 3 *) 2 *3))) / 

(4({-l+E J “ a ) 2 -4E? h3<, h3V)) + 

(mg (-3+K2) - a 2 (1+X2)) . 

8(-y + n*.) 

042= (E?’ 211 ^ (-a (-1+ E! 41130 + 4ET 2h3 “h3a) <l + x 3 ) + 

n fe (l+EJ 4h3a -2I^“ a (l+4h3V) - (-1 +E? KJo ) 2 x 3 )))/ 
(4((-l+E 2Ma ) 2 -4E?““h3 2 a 2 )) + 

E? 0 * 2 (m? (-3 + X2) - a 2 (1+ *2) ) . 

8 <-Y +Hte) 

C43, (E* 2 "* (-a (-l*S 4Ua + 4E 2Mo h3a) (l+x 3 ) + 

(1+ E 42 **" - 2E? h3a (1+ 4h3 2 a 2 ) - (-1+ E?“ a ) 2 x 3 ) ) ) / 
(4 ((-1+ E? h3a ) 2 -4E? h3a b3V}) + 

(mj| {-3+ X2) - a 2 (1+ X2) ) t 

8(-y+iife) 

044 = (&*"* (-a (- 1+ El 4h3 “ + 4 E? 11 * 0 h3a) (l+x 3 ) + 

*, (!♦ - 2E?“ a (1+ 4h3 2 a 2 ) - (-1+ E? M “) 2 x 3 ) ) ) / 

(4((-l+E?“ a ) 2 -4E 2 “ a h3?a 2 )) + 

E? 0 * 1 (ir^ ( -3 + x 2 ) - a 2 { 1+ x2) ) . 

8 (-r +in») 
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